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\ Abstract. In order to investigate the formation of O I 7771-5 and [O i] 6300/6363 

i ! 

, lines, extensive non-LTE calculations for neutral atomic oxygen were carried out 



for wide ranges of model atmosphere parameters, which are applicable to early-K 
through late-F halo/disk stars of various evolutionary stages. 

The formation of the triplet O I lines was found to be well described by the clas- 
sical two-level-atom scattering model, and the non-LTE correction is practically 
q . determined by the parameters of the line-transition itself without any significant 

■ relevance to the details of the oxygen atomic model. This simplifies the problem in 

CZ2 . 

the sense that the non-LTE abundance correction is essentially determined only 

by the line-strength (W\), if the atmospheric parameters of T e a, \ogg, and £ are 
■ *~ i ■ 

given, without any explicit dependence of the metallicity; thus allowing a use- 

;_h ' 

ful analytical formula with tabulated numerical coefficients. On the other hand, 
our calculations lead to the robust conclusion that LTE is totally valid for the 
forbidden [O i] lines. 

An extensive reanalysis of published equivalent-width data of O I 7771-5 and 
[O i] 6300/6363 taken from various literature resulted in the conclusion that, 
while a reasonable consistency of O I and [O i] abundances was observed for disk 
stars ( — 1 [Fe/H] 0), the existence of a systematic abundance discrepancy 
was confirmed between O I and [O i] lines in conspicuously metal-poor halo stars 
(—3 ^ [Fe/H] —1) without being removed by our non-LTE corrections, i.e., the 
former being larger by ~ 0.3 dex at —3 ^ [Fe/H] 5^ —2. 

An inspection of the parameter-dependence of this discordance indicates that 
the extent of the discrepancy tends to be comparatively lessened for higher 
T e ff / log g stars, suggesting the preference of dwarf (or subgiant) stars for studying 
the oxygen abundances of metal-poor stars. 
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1. Introduction 

Oxygen abundance determination in F-K stars is one of the most controversial and 
interesting topics in stellar spectroscopy, especially concerning the behavior of the [O /Fe] 
ratio in metal-poor halo stars (see, e.g., Carretta et al. 2000; King 2000; Israelian et al. 
2001; Nissen et al. 2002; and the references therein). 

Namely, we are faced with a problem currently in hot controversy on the observational 
side, which stems from the fact that each of the different groups of spectral lines, (a) 
high-excitation permitted O I lines (e.g., 7773 triplet), (b) low-excitation forbidden [O i] 
lines (6300/6363 lines), and (c) OH lines at UV (~ 3100-3200 A), do not necessarily 
yield reasonably consistent oxygen abundances in metal-poor stars. 

More precisely, a distinct discrepancy is occasionally reported between the abundances 
derived from O I and [O i] lines in the sense that the former tend to yield systematically 
higher oxygen abundances than the latter in metal-poor stars (see, e.g., Cavallo et al. 

1997) , in contrast to the relation between O I and OH lines which give more or less 
consistent abundances (e.g., Boesgaard et al. 1999). As a result, different trends of O-to- 
Fe ratio, an important tracer of early nucleosynthesis in our Galaxy, have naturally been 
suggested: (1) plateau-like [O/Fe] from [O i] lines, and (2) ever-rising [O/Fe] from O I or 
OH lines. 

Admittedly, what has been described above is nothing but a rough summary and the 
actual situation is more complicated than such a simple dichotomic picture. For example, 
the O I vs. [O i] discordance appears to be observed only in metal-poor halo stars, since 
such a systematic disagreement is absent in population I G-K giants (cf. Takeda et al. 

1998) . Moreover, according to Carretta et al.'s (2000) recent investigation of very metal- 
poor stars (mainly giants), a "nearly flat" [O/Fe] trend was concluded consistently from 
both of the O I and [O i] lines. Very recently, on the other hand, Nissen et al. (2002) also 
obtained nearly consistent O I and [O i] abundances of metal-poor dwarfs and subgiants, 
but at a "quasi- linearly increasing" trend in contrast to Carretta et al. (2000), though 



* Tables 2, 5, and 7 are available only in electronic form, which are temporarily at the anony- 
mous ftp site of ftp://133.ll.160.242/Users/takeda/aams2727/ but will be replaced (sub- 
mitted) to CDS when the paper is to be published. 
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an application of the 3D hydrodynamical model atmospheres 1 (taking into account the 
atmospheric granular motions) tends to deteriorates the agreement. 

This confusion makes us wonder which on earth is the correct and reliable indicator of 
the oxygen abundance. Generally speaking, a modest view currently preferred by many 
people is that [O i] lines are less problematic and more creditable due to their weak 
T e ff-sensitivity as well as the general belief of the validity for the assumption of LTE, 
while O I and OH lines are less reliable because of the uncertainties involved with the 
non-LTE correction (O i) or with the 3D effect of atmospheric inhomogcncity (OH). 

Recently, however, we cast doubt on this conservative picture and suggested the 
possibility of erroneously underestimated abundances from [O i] lines (i.e., the problem 
may be due to the formation of forbidden lines rather than that of permitted lines) 
based on our analysis of seven metal-poor stars (Takeda et al. 2000), which indicated the 
tendency of above-mentioned O I vs. [O i] discrepancy (though a small number of sample 
stars prevented us from making any convincing argument). 

Motivated by that work suggesting the necessity of a follow-up study to check if our 
hypothesis is reasonable or not, we decided to revisit this O I vs. [O i] problem while 
paying special attention to the topic of the non-LTE correction to be applied, one of 
the controversial factors in oxygen abundance determination. Toward this aim, we car- 
ried out new extensive non-LTE calculations, based on which the numerous published 
equivalent-width data were reanalyzed to derive the non-LTE abundances, in order to dis- 
cuss whether or not the abundance discrepancy between O I and [O i] really exists, which 

1 This effect is another controversial factor of importance in the oxygen abundance problem, 
and the results from the classical ID model atmospheres neglecting this effect (e.g., the tem- 
perature of the upper layer is considerably lowered in the 3D atmospheric model) may suffer 
appreciable corrections. Though theoretical investigations in this field are still in progress, the 
nature of the 3D abundance corrections ( A3D log e) to be applied to the ID solutions for each 
oxygen line may be roughly summarized as follows (e.g., Kiselman & Nordlund 1995; Asplund 
2001; Allende Prieto et al. 2001; Asplund & Garcia Perez 2001; Nissen et al. 2002): 

— This 3D effect is practically negligible for the deeply-forming high-excitation O I lines (cf. 
Fig. 10 in Kiselman & Nordlund 1995; Asplund 2001). 

— Similarly, its effect on the forbidden [O 1] line is generally minor (IA3D log e| ^ 0.1 dex; the 
sign of the correction is negative), but may be increased to an appreciable amount up to ~ 0.2 
dex in lower gravity /met allicity atmospheres (cf. Table 6 in Nissen et al. 2002). 

— The 3D effect is most important on the UV OH lines, and conspicuously large corrections 
amounting to |A3Dloge| ~ 0.6 dex or even more (at [Fe/H] = —3.0; Asologe being negative) 
are suggested (cf. Table 2 in Asplund & Garcia Perez 2001). 

Since our main concern in this paper is to investigate the non-LTE effect on the abundances 
derived from the O I and [O 1] lines, our use of the classical ID models would not cause any 
serious problems (actually, taking account of the 3D correction does not essentially affect our 
final conclusion; cf. the footnote 11 in Sect. 4.5.1). 



4 Y. Takcda: Oxygen line formation in late-F through early-K disk/halo stars 

we hope to shed some light on the oxygen problem of metal-deficient stars mentioned 
above. 

This paper is organized as follows: 

— First, extensive tables of non-LTE corrections for the O I 7771-5 triplet lines over 
a wide range of atmospheric parameters are presented, so that they may be applied to 
late-F through early-K dwarfs/giants of various metallicities. 

— Second, in order for application of such tables to actual analysis of observational data, 
we try to find out a useful analytical formula for evaluating non-LTE corrections of these 
O I triplet lines for any given atmospheric parameters (T e ff, log <?, £) and an observed 
equivalent width (W\). 

— Third, we investigate the possibility of a non-LTE effect for the low-excitation forbid- 
den [O i] 6300/6363 lines, for which few such attempt has ever been made. 

Fourth, equivalent-width data of O I 7771-5 lines as well as [O i] 6300/6363 lines taken 
from various literature are reanalyzed within the framework of the consistent system of 
analysis along with application of our non-LTE correction formula, in order to discuss 
the characteristics of [O/Fe] vs. [Fe/H] relation derived from permitted O I lines and 
forbidden [O i] lines separately. 

Also, as a supplementary subject related to this study, the influence of the treatment of 
H I collision on the non-LTE corrections of O I 7771-5 lines was quantitatively estimated, 
and the consistency of the present choice of this parameter was checked by comparing 
the abundances of nearby solar-type stars derived from O I 7771-5, O I 6158, and [O i] 
6300 lines. In addition, we also tried to qualitatively explain the line formation of O I 
7771-5 in terms of the simple classical two-level-atom model, which gives us an insight 
to understanding the behavior of the non-LTE correction for these triplet lines. These 
are separately described in Appendices A and B, respectively. 

2. Statistical equilibrium calculations 

2.1. Atomic model 

The procedures of our non-LTE calculations for neutral atomic oxygen are almost the 
same as those described in Takeda (1992), which should be consulted for details. We only 
mention here that the calculation of collision rates in rate equations basically follows the 
recipe described in Takeda (1991). Especially, the effect of H I collisions, an important 
key factor in the non-LTE calculation of F-K stars, was treated according to Steenbock 
& Holweger's (1984) classical formula without any correction as has been done in our 
oxygen-related works so far, which we believe to be a reasonable choice as far as our 
calculations are concerned (cf. Takeda 1995, and Appendix A). 

There is, however, an important change concerning the atomic model as described 
below. The model oxygen atom adopted in Takeda (1992) was constructed mainly based 
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on the atomic data (energy level data and gf values) compiled by Kurucz & Peytremann 
(1972), with ~ 30 additional radiative transitions which were missing in their compilation 
[e.g., resonance transitions whose transition probabilities were taken from Morton (1991)] . 
This old model did not include the 2s 2 2p 4 3 P-2p 4 1 D transition (corresponding to 
[O i] 6300/6363 lines) explicitly as the radiative transition (i.e., it was treated as if 
being completely radiatively forbidden in spite of its actually weak connection), simply 
because it was not contained in Kurucz & Peytremann (1972). We also noticed that 
several important UV radiative transitions originating upward from the 2p 1 D term 
(i.e., upper term of [O i] 6300/6363) were also missing in that old model atom for the 
same reason. 

We, therefore, newly reconstructed the model oxygen atom exclusively based on the 
data of Kurucz & Bell (1995), with additional transition probability data for 2p 4 1 D- 
3s 3 S° and 2p 4 1 D-3s 5 S° transitions taken from Vienna Atomic Line Database (VALD; 
Kupka ct al. 1999). Regarding the collisional rate Ci_2 for the [O i] transition (2s 2 2p 4 
3 P-2p 4 1 D) mentioned above, the cross-section given in Osterbrock (1974) were used for 
evaluating the electron collision rates (C 1 _ 2 ), and the H I collision rates (C{ 1 _ 2 ) were 
estimated by scaling Cf_ 2 (cf. Takeda 1991). The resulting atomic model comprises 87 
terms (including terms up to 3p" 1 S, 6d' 3 P, lOd 5 D for the singlet, triplet, and quintet 
system, respectively) and 277 radiative transitions, which is similar to the old one in 
terms of the complexity but must be significantly improved, especially for studying the 
formation of [O i] lines (i.e., more realistic treatment of the 1-2 interaction and the 
inclusion of as many UV transitions connecting to the upper term 2 as possible). 

Given the changes explained above, the reader can find the description on the adopted 
rates in the statistical-equilibrium calculation in Sect. 2.1 of Takeda (1992). Some addi- 
tional details (which were not explicitly described therein) are given below: 

— The population of neutral hydrogen (necessary for computing the rates of H I colli- 
sion or the charge-exchange reaction) was taken from the model atmospheres (i.e., LTE 
ionization equilibrium). 

— For the upper terms other than the lowest 8 terms mentioned in Sect. 2.1 of Takeda 
(1992), the hydrogenic approximation was used for the photoionization cross section. 

— The photoionizing radiation field was computed based on Kurucz's (1993a) ATLAS9 
model atmospheres while incorporating the line opacity with the help of Kurucz's (1993b) 
opacity distribution function. 

— The collisional rates for those transitions, which are not explicitly mentioned in Sect. 
2.1 of Takeda (1992), were evaluated according to the procedure described in Sect. 3.1.3 
of Takeda (1991). 
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2.2. Model atmospheres 

We planned to make our calculations applicable to stars from near-solar metallicity (pop- 
ulation I) down to very low metallicity (extreme population II) at late-F through early-K 
spectral types in various evolutionary stages (i.e., dwarfs, subgiants, giants, and super- 
giants). We, therefore, decided to carry out non-LTE calculations on the extensive grid 
of one hundred (5x5x4) model atmospheres resulting from combinations of five T c s 
values (4500, 5000, 5500, 6000, 6500 K), five logg (cm s~ 2 ) values (1.0, 2.0, 3.0, 4.0, 5.0), 
and four metallicities (represented by [Fe/H]) (0.0, —1.0, —2.0, —3.0). As for the stellar 
model atmospheres, we adopted Kurucz's (1993a) ATLAS9 models 2 corresponding to the 
microturbulent velocity (£) of 2 kms -1 . 

2.3. Oxygen abundance and microturbulence 

Regarding the oxygen abundance used as an input value in non-LTE calculations, we 
assumed loge£ put = 8.93 + [Fe/H] +[0/Fc], where [O/Fe] = 0.0 for the solar metallicity 
models ([Fe/H] = 0) and [O/Fe] = +0.5 for the other metal-deficient models ([Fe/H] 
= —1, —2, and —3). Namely, the solar oxygen abundance of 8.93 taken from Anders 
& Grevesse (1989), which is also the canonical value used in the ATLAS9 models, was 
adopted for the normal-metal models, while the metallicity-scaled oxygen abundance plus 
0.5 dex (allowing for the characteristics of a element) was assigned to the metal-poor 
models. 

The microturbulent velocity (appearing in the line-opacity calculations along with 
the abundance) was assumed to be 2 kms" 1 , to make it consistent with the model 
atmosphere. 



2 It has been occasionally argued that the new treatment of convection in these ATLAS9 
models may not be adequate (e.g., Castelli et al. 1997) in the sense that the newly incorporated 
effect of convective overshooting had better be rather switched off [like the old Kurucz's (1979) 
ATLAS6 models]. In order to investigate this problem, test calculations were carried out by 
using the models without convective overshooting (generated by the ATLAS9 program while 
setting OVERWT = 0) and the abundance variations between the overshooting-on and -off 
cases were examined for a wide range of atmospheric parameters (5000 K < T e g < 6500K, 
3 < log5 < 5, and -3 < [Fe/H] < 0). It turned out that the differences of log e( ATLAS9 model 
with overshooting) — log e(no-overshooting model) are mostly positive for both cases of O I and 
[O i], and quantitatively insignificant; i.e., ^ 0.05 dex for [O i] and Ss 0.1 dex for O I (which are 
interpreted as reflecting the change in the temperature gradient). Considering the same sign of 
these variations, we may consider that this effect can be neglected in the present O I vs. [O i] 
problem. 
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3. Non-LTE corrections 

3.1. Evaluation procedures 

Here, we describe how the non-LTE abundance corrections are evaluated after the 
departure coefficients, b = Hnlte / ^lte > have been established by statistical equilibrium 
calculations. 

For an appropriately assigned oxygen abundance (A a ) and microturbulence (£ a ), we 
first calculated the non- LTE equivalent width (W NLTE ) of the line by using the computed 
non-LTE departure coefficients (6) for each model atmosphere. Next, the LTE (A L ) and 
NLTE (A N ) abundances were computed from this Fy NLTE while regarding it as if being a 
given observed equivalent width. We can then obtain the non-LTE abundance correction, 
A, which is defined in terms of these two abundances as A = A N — A L . 

Strictly speaking, the departure coefficients [b(r)] for a model atmosphere correspond 
to the oxygen abundance and the microturbulence of log£Q PUt and 2 kms -1 adopted 
in the calculations (cf. Sect. 2.3). Nevertheless, considering the fact that the departure 
coefficients (i.e., ratios of NLTE to LTE number populations) are (unlike the population 
itself) not much sensitive to small changes in atmospheric parameters, we applied such 
computed b values also to evaluating A for slightly different A & and £ a from those fidu- 
cial values assumed in the statistical equilibrium calculations. 3 Hence, we evaluated A 
for three A & values (log£Q PUt and ±0.3 dex perturbation) as well as three £ values (2 
kms -1 and ±1 km s _1 perturbation) for a model atmosphere using the same departure 
coefficients. 

3.2. Results of computed corrections 

The computation of A values was carried out for five lines: three O I lines of 7771-5 
triplet and two [O i] lines of 6300/6363 doublet. The WIDTH9 program, which was 
originally written by R. L. Kurucz but modified by the author in various respects (e.g., 
incorporation of the non-LTE departure coefficients into the line opacity and the line 
source function), was used for calculating the equivalent width for a given abundance, or 
inversely evaluating the abundance for an assigned equivalent width. The adopted line 
data are summarized in Table 1. 

3 The error caused by this practical treatment was checked based on exact calculations for the 
representative case of T eS = 6000 K, log 3 = 3.0, and [O/Fe] = [Fe/H] = 0.0, where the non-LTE 
correction is rather large (A ~ —0.5 dex). It was confirmed that the [correct treatment — prac- 
tical treatment] difference in evaluating the non-LTE correction is quantitatively insignificant; 
i.e., 0.01 dex [the use of b (2 km s -1 ) for computing A(2 ± 1) km s _1 ] and 0.04 dex [the use 
of 6(logeQ PUt ) for computing A(logeQ PUt ± 0.3)]. 
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The full results of the computed A values (£ = 1, 2, and 3 kms" 1 cases for each of 
the three triplet lines at 7771.94, 7774.17, and 7775.39 A) are given in Table 2. 4 

As will be described in Sect. 4.3, the non-LTE corrections for the forbidden [O i] 
6300/6363 lines turned out to be essentially equal to zero, suggesting the firm validity of 
LTE. 

3.3. Analytical formula for practical application 

It is necessary for us to evaluate the non-LTE corrections corresponding to actual 
equivalent- width data for stars having a wide variety of atmospheric parameters (T e s, 
log 5, £) and metallicity ([Fe/H]), in order to compare them with those derived by others, 
or to apply them to reanalysis of published data. However, since the raw tabulated values 
given for each of the individual models (e.g., Table 2) are not necessarily convenient for 
such practical applications, it is desirable to develop a useful formula for this purpose. 

This motivation reminded us of the simple unique correlation between W\(0 I 
7774.17) and its NLTE correction in metal-poor solar-type dwarfs (5500 K < 
T cS < 6500K, 3.5 < log 5 < 4.5, and -2 < [Fe/H] < 0) discussed in Takeda (1994; 
cf. Fig. 10 therein). 5 If such a relation really holds also in the present case irrespective 
of the atmospheric parameters (much wider range than that in Takeda 1994), it would 
greatly simplify our task. However, according to the A vs. FT 1 ^ 1 "™ relation for O 1 7774.17 
depicted in Fig. 1 (based on the data in Table 2), the situation is not such simple; i.e., 
A is not only a function of FT NLTE but also more or less dependent 6 on T e s and logg 
(and £). Nevertheless, a close examination revealed that for given T e g, log g, and £, the 
extent of the non-LTE correction (A) is a nearly monotonic function of the equivalent 
width (W\; measured in mA), and that 

A = al0 6w \ (1) 

where a and b are appropriately chosen coefficients as functions of (T e g, log 5, and £), 
is a fairly good approximation for an appropriately chosen set of (a, b), as far as the line 
is not too strong (i.e., W\ ^ 100mA). It should be noted that this relation does not 
explicitly contain the metallicity, though stars with less metals generally show smaller 
extent of |A| values due to their weaker line-strengths (W\). 

4 Available only in electronic form (cf. the footnote to the title). Similarly, the non-LTE 
corrections for another line, O I 6158, were calculated for the [Fe/H] = and —1 models for the 

discussion in Appendix A, and are presented as (electronic) Table 7. 

5 This tendency can be interpreted as being due to the close connection between the extent 

of the dilution of the line source function (determining the extent of the non-LTE effect) and 

the strength of these triplet lines (cf. Sect. 4.2 in Takeda 1994, or Appendix B). 

6 This dependence on T c g and logp may be qualitatively explained by the simple two-level- 
atom model, as shown in Appendix B. 
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We therefore determined the best-fit coefficients (a, b) based on the computed A 
values in Table 2 for each combination of (T e ffj log<7, £); the resulting coefficients are 
summarized in Tables 3 for each of the O I triplet lines at 7771.94, 7774.17, and 7775.39 
A, respectively. The analytical approximations presented by Eq. (1) with the coefficients 
given in Table 3 are also depicted in Fig. 1 (the solid lines). 

Practically, we first evaluate (a, b) by interpolating this table in terms of T e fj, logg 
and £ of a star in consideration, which are then applied to Eq. (1) with the observed 
equivalent width to obtain the non-LTE correction. 

In order to demonstrate the applicability of this formula, we show the error of Eq. 
(1) (defined as 5 mor = A formula - A actual ) as a function of W\ for the case of O I 7774.17 
in Fig. 2 (compare it with Fig. 1), where we can see that an accuracy within a few 
hundredths dex is accomplished as far as W\ Si 100 mA. 

4. Discussion 

4.1. NLTE correction for O I triplet: comparison with others 

It is interesting to compare the non-LTE corrections for the O I 7771-5 reported in 
several published studies with those estimated following the procedure described in Sect. 
3.3 while using the same observational equivalent widths along with the same T c ff , logy 
and £ as given in the literature. If W\ values for two or three lines of the triplet are 
available, we computed the corrections for each of the lines and their average was used 
for the comparison. 

4.1.1. Takeda et al. (1998, 2000) 

Comparisons with the results of our own previous work, Takeda et al. (1998, 2000), are 
shown in Figs. 3a and b. 

The non-LTE corrections for O I 7771.94 given in Table 2 of Takeda et al.'s (1998) 
population I G-K giants analysis are in good agreement with those derived from our ap- 
proximate formula [Eq. (1)], in spite of the fact that extrapolating Table 3 was necessary 
for stars with (4250 K <,)T cS < 4500 K. 

We can also see that the A values of O I 7771-5 for very metal-poor giants computed 
by Takeda et al. (2000) are reasonably consistent with those estimated by the procedure 
proposed in this study. Note, however, that a somewhat large discrepancy is observed 
in the £ = 5.1 kms -1 case for HD 184266, where two different £ values were suggested 
depending on the lines (O I triplet lines and Fe I lines) used for its determination. This 
indicates that we should be careful at extrapolating Table 3 to considerably outside of 
the range they cover (e.g., 1 kms _1 < £ < 3 kms -1 ). 
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Generally speaking, the results (for O I lines) of our new calculations reported in 
this paper (cf. Table 2) are essentially equivalent to our previous studies, in spite of the 
use of updated new atomic model in connection with the [O i] transition (cf. Sect. 2.1). 
This can be verified also by comparing Table 2 in this paper with Table 7 of Takeda 
(1994), both of which are practically equivalent. Taking this fact into consideration, we 
can again confirm from the coincidence shown in Figs. 3a and b that Eq. (1) with Table 
3 is sufficiently accurate and useful for practical applications. 

4.1.2. Tomkin et al. (1992) 

In Tomkin et al.'s (1992) extensive analysis of C and O abundances in (from mildly to 
very) metal-deficient dwarfs, they performed non-LTE calculations by using a 15 levels/22 
lines atomic model while adjusting the H I collision rate by the requirement that the solar 
O I 7771-5 lines yield the abundance of 8.92 with the solar model atmosphere of Holweger 
& Mullcr (1974). As shown in Fig. 4a, the extents of their |A| values (< 0.1 dcx) tend 
to be systematically smaller than ours. 

4.1.3. Mishenina et al. (2000) 

Comparing our A values with the recent non-LTE calculation of Mishenina et al. 
(2000), who used an atomic model comprising 75 levels/46 transitions and Steenbock 
& Holweger's (1984) H I collision formula with a reduction factor of 1/3, we can see that 
both arc in reasonably good agreement in spite of a mild difference in the treatment of 
H I collision rates (cf. Fig. 4b). 

4.1.4. Nissen et al. (2002) 

Nissen et al. (2002) recently carried out a careful abundance study based on their high- 
quality data of O I and [O i] lines (along with OH lines) , and obtained almost consistent 
results between [O i] (LTE) and O I (non-LTE) at a "quasi-linearly increasing" trend 
of [O/Fe] within the framework of the ID atmosphere analysis (though an application 
of the 3D correction to the [O i] abundance causes some discrepancy). Their O I non- 
LTE calculation is based on an atomic model comprising 22 levels/43 transitions, while 
completely neglecting the H I collisions. We can see from Fig. 4c that the extent of their 
|A| values tends to be larger than ours by 5s 0.15 dex, which may be understood as being 
due to the different treatment of H I collisions, since neglecting this effect increases |A| 
by this amount (cf. Appendix A). Even with this moderate difference, the reanalysis of 
their W\(0 i) and Wa([0 i]) data based on our A values nearly reproduces the results 
they obtained (i.e., the similar tendency for O I and [O i]), as shown in Fig. 8a. 
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4.1.5. Carretta et al. (2000) 

Carretta et al. (2000) also obtained a nearly consistent solution of [O/Fe] between O I 
and [O i] in the metal-poor regime, but at a "nearly flat" trend in contrast to the results 
of Nissen et al. (2002) mentioned above. Namely they reanalyzed the published W\(0 I 
7771-5) data of Tomkin et al. (1992) and Edvardsson et al. (1993) and showed that their 
(NLTE-corrected) O I abundances are consistent with those of [O i] resulting from their 
reanalysis of W x ([0 i] 6300/6363) data of Snedcn et al. (1991) and Kraft et al. (1992). 
Their corrections are based on Gratton et al.'s (1999) NLTE calculations, where they 
adopted a model oxygen atom of 13 levels/28 transitions and a mild enhancement factor 
of <~ 3 (determined by the requirement of O I 7771-5/0 I 6158 consistency in RR Lyr 
variables) applied to classical H I collision rates. 

In Fig. 5a, their NLTE corrections are compared with our corresponding A val- 
ues, which were evaluated from the W\(0 I 7771-5) data of Tomkin et al. (1992) and 
Edvardsson et al. (1993) while adopting the T e s and log g values used by Carretta et 
al. (2000) to make the comparison as consistent as possible. An interesting tendency 
is observed in this figure. That is, while the extents of their non-LTE corrections for 
Edvardsson et al.'s (1993) F-G disk dwarfs are systematically smaller than ours, this 
trend disappears and a good agreement is seen in the case of Tomkin et al.'s (1992) 
considerably metal-poor halo dwarf stars (their |A| values are even larger for two stars). 

A closer inspection revealed that this is nothing but the dependence on the metallicity, 
as shown in Fig. 5b. Namely, their corrections (in terms of the absolute |A| values) are 
smaller than ours for mildly metal-deficient disk stars (—1 [Fe/H] ^ 0), nearly the 
same for typical halo stars (—2 ^ [Fc/H] ^ — 1), and eventually become larger than our 
values for very metal-poor halo stars (—3 ^ [Fc/H] ^ —2). Evidently, this must be the 
reason for their "quasi-flat" trend of [O/Fc] obtained for the O I lines, since the extent 
of their (negative) non-LTE correction becomes progressively larger than ours (in the 
relative sense) as the metallicity decreases, thus suppressing the gradient of the [O/Fe] 
vs. [Fc/H] relation. 

4.2. Does the metallicity affect the formation of O I lines? 

In Sect. 4.1.5 we found a metallicity-dependent systematic discrepancy between Carretta 
et al.'s (2000) non-LTE corrections and ours (i.e., theirs are systematically overcorrected 
toward a lowered metallicity compared to ours), which must be a very important problem 
because of its large impact on the [O/Fe] vs. [Fe/H] relation derived from the O I triplet 
lines. Hence, a detailed discussion on this matter from the viewpoint of the line-formation 
mechanism may be due. 

Actually, such a [Fc/H]-dependent tendency of their correction is clearly observed 
by examining Table 10 of Gratton et al. (1999). We plotted their A vs. W A (7771.94) 
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relations computed for the log g — 3.0 models of different mctallicity in Figs. 6a and b, 
each corresponding to T e g = 5000 K and 6000 K, where our computed results (Table 
2) along with the corresponding analytical approximation curves [i.e., Eq. (1) with the 
coefficients taken from Table 3] are ovcrplottcd. ft can be seen from these figures that 
the strong mctallicity-dependence mentioned above is evidently observed in their results 
(especially for the case of T c s = 6000 K). However, such a tendency (i.e., A is considerably 
dependent on the metallicity at the same W\) is not seen in our case, where A is nearly 
a monotonic function of W\ as a first approximation (see in Figs. 6a and b that the small 
symbols of different type continuously follow a common curve with a change in W\). 

In order to demonstrate this more clearly, additional test non-LTE calculations were 
performed for the T e g = 6000 K and logg = 3.0 models of different metallicities 
([X] = 0,-1,-2,-3) where the assigned oxygen abundances were carefully adjusted 
so as to yield the same (non-LTE) O I 7771.94 equivalent widths for the two metallicity 
pairs, W x = 152 mA (for [X] = and -1) and W\ = 20 mA (for [X] = -2 and -3) 
[which were chosen to keep consistency with Gratton et al.'s (1999) calculation], and 
computed the non-LTE correction for each case. The results are given in Table 4, where 
the corresponding corrections taken from Gratton ct al.'s (1999) Table 10 are also quoted 
for a comparison. It is clear from this table that such a marked metallicity-dcpcndence 
observed in Gratton ct al.'s (1999) corrections (column 12) is absent in our A values 
(column 11). 

The reason for this characteristics in our case (i.e., near-equality of A at the same 
W\ but different metallicities) is manifestly explained by Fig. 7, where the behaviors of 
the source function, the NLTE-to-LTE line-opacity ratio, and the Planck function for the 
four test cases of Table 4 are depicted. We should note the following facts which can be 
read from this figure: 

(1) The line opacity for this O I line is nearly in LTE (^ LTE /# TE ~ 1) in the line- 
forming region. 7 Though an appreciable departure from LTE (i.e., overpopulation) is 
observed at the shallower layer especially for the metal-poor atmospheres, it does not 

7 This is closely related to the fact that oxygen atoms are dominantly in the neutral stage 
in the line-forming region within the framework of our classical model atmosphere, where LTE 
essentially holds in the O i/O n ionization equilibrium due to the efficient charge-exchange 
reaction with hydrogen atoms (whose population/ionization is computed on the assumption of 
LTE). [For example, even in the high-temperature, low-gravity, very metal-poor atmosphere of 
(T e ff = 6500 K, log g — 1, [Fe/H] = —3) where the ionization is considered to be most important 
among our studied cases, ~ 99% of oxygen is neutral at T5000 ~ 0.1.] Accordingly, the lower-level 
population of the O I triplet is significantly affected by the population of the ground level where 
most of the O I atoms populate and LTE almost prevails. Note that the photoionizing radiation 
field (which is conspicuously dependent on the metallicity via the line/continuum opacities) does 
not play any important role for the O I population. 
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have any significant affection upon the emergent radiation. [See Sect. 4.2 in Takeda (1994) 
for the explanations on these characteristics of ^ lte /^q TE .] 

- (2) As can be seen from Fig. 7, the line source function (Sl) of this O I 7771.94 
line resulting from our detailed non-LTE calculations (thin solid lines) is essentially the 
same as the two- level-atom- type line source function (open circles) . Namely, the solution 
obtained by solving the classical scattering problem with the simple form of ,g£ wo_level 5 

£two-lcvcl = g + e '£)/(! + £ ') ; ( 2 ) 

while using the realistic values of e'(r), B(t), and ij(t) computed from the model atmo- 
spheres [where J (frequency-averaged mean radiation field), e' (collisional-to-radiative 
ratio of the deexcitation rate), B (Planck function), and r/(r) (the line-to-continuum 
opacity ratio necessary to compute the line optical depth and to combine the line-and- 
continuum source function into the total source function) are the quantities related to 
the transition between the two levels; cf. Eq. (11-6) in Mihalas (1978) or Appendix B 
in this paper], yields essentially the same result as S^" 10 , obtained from the solutions of 
detailed statistical-equilibrium calculations, 

Sl™ c = {2h^/c 2 )[{n l g j )/(n, J g l ) - l]" 1 , (3) 

where rn/rij and gi/gj are the number population and the statistical weight of 
lower/upper level, respectively [cf. Eq. (4-14) in Mihalas (1978)]. This evidently sug- 
gests that the formation of the I triplet is essentially determined by the quantities 
directly related to this transition itself without any significant relevance to the details of 
the adopted atomic model. 8 

- (3) In this the well-known characteristics of such a classical "two-level-atom" 
problem (e.g., Hummer 1968; Appendix B), the extent in the dilution of the line source 
function (S^/B) is controlled by the line-to-continuum opacity ratio n (= Io/k). And 
since both eventually determine the strength of the line, it is quite natural that the non- 
LTE effect (i.e., the extent of the S*l dilution) is closely coupled with the line strength 
(Wx). 

Then, there is almost no room for the metallicity directly playing any significant role 
in the non-LTE effect. Admittedly, there are some possibilities of metallicity-dependence 
even within the framework of such a two-level-type line-formation. However, any of 



8 The line source function of any transition [given by Eq. (3)] can be generally expressed (by 
rewriting the rate equations) as S^ uc = (J + e' B + ci)/(l + e' + C2), where ci and C2 are the 
interlocking terms representing the cumulative influences of other levels and transitions [cf. Eq. 
(12-4) in Mihalas (1978)]. Then, the consequence of Sf uo ~ S^°~ level naturally implies that 
the complicated effects of other levels and transitions (represented by ci and C2) are almost 
negligible; i.e., these O I 7771-5 lines form as if the oxygen atom were made of only the two 
levels relevant to this transition. 



14 Y. Takeda: Oxygen line formation in late-F through early-K disk/halo stars 

them are unlikely in the present case. For example, appreciable [Fc/H]-dcpcndences of T 
(metallicity-related surface cooling) or of ^ lte /^q TE [see (1) above] existing near to the 
surface (cf. Fig. 7) do not matter, because they occur far above the line-forming region. 
Also, the collisional rate between two levels (affecting the photon destruction probability 
e') is hardly dependent on the mctallicity (or electron density) in our case because of 
the dominance of the H I collisions 9 [though an extreme case of pure electron collisions 
(neglecting H I collisions) would lead to a metallicity-dcpendence of the e' values (and 
eventually of the non-LTE effect) owing to the role of mctallicity as the electron donor; 
e.g., Kisclman (1991)]. 

A supplementary discussion based on the two-level-atom along with the simple Milne- 
Eddington model is given in Appendix B, where the T e ff/logg-dependence of the A vs. 
W\ relation given in Fig. 1 (along with its [Fe/H]-independence) is qualitatively explained 
more in detail. 

Accordingly, the results of the Padova group, that the extent of the non-LTE correc- 
tion increases with a decrease in the metallicity (for a given W\), is hard to understand at 
least from the viewpoint of what we found from our calculations. Namely, if their results 
are to be seriously taken, something in their calculations must be greatly different from 
ours (e.g., breakdown of the two-level-atom nature and the dominant contributions from 
other levels and transitions). Curiously, however, an inspection of Gratton et al.'s (1999) 
Fig. 6 suggests that their non-LTE departure coefficients behave themselves consistently 
with our computational results (see, e.g., Fig. 9 of Takeda 1994); that is, for the case 
of T cff = 6000 K and \ogg = 4.5, the departure from LTE is larger for the [Fc/H] = 
model than for the [Fe/H] = —2 model. In other words, their non-LTE calculations 
appears quite reasonable in this respect. Then, how could their calculations lead to the 
tendency of enhanced | A| with a lowered [Fe/H] ? Anyway, it appears that the non-LTE 
corrections given in Gratton et al.'s (1999) Table 10 and those used by Carretta et al. 
(2000) are questionable, and thus their consequences had better be viewed with caution. 
This criticism also applies to any other results obtained by using their corrections; e.g., 
the "flat-like" [O/Fe] tendency derived by Primas et al. (2001). 



9 For example, the values of Ch/C (at tbooo ~ 0.1) for the O I 7771-5 transition, computed 
from Steenbock &Holweger's (1984) formula (for Ch) and Van Regemorter's (1962) formula (for 
Co) (as adopted in this study) for each of the representative models with parameters of (T e fr, 
logp, [Fe/H]) are as follows: ~30/~1000 for (4500 K, 4.0, 0/-2), ~8/~500 for (4500 K, 1.0, 
0/-2), ~7/~100 for (5500 K, 4.0, 0/-2), ~4/~10 for (5500 K, 1.0, 0/-2), ~2/~6 for (6500 K, 
4.0, 0/-2), and ~0.3/~0.5 for (6500 K, 1.0, 0/-2). Hence, except for the high-temperature and 
low-gravity model (which is almost irrelevant in the present study of metal-poor stars), Ch is 
generally predominant over C e . 
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4.3. Validity of LTE for the [0 I] 6300/6363 lines 

As mentioned in Sect. 1, the search for the possibility of the non-LTE effect in the [O i] 
lines, based on the new (more realistic) O I atom, was actually one of the main aims in 
this paper. 

According to the present calculation, however, the departure coefficients of the lower 
term (2s 2 2p 4 3 P; hereinafter referred to as "term 1") and the upper term (2p 4 1 D; 
hereinafter "term 2") for the [O i] 6300/6363 lines turned out to be essentially equal to 
unity over almost all atmospheric layers, which guarantees the validity of the assumption 
of LTE for these forbidden [O i] lines. 

The reason for this is attributed to the fact that (1) the population of the lower ground 
term is thermalized relative to the continuum (i.e., accomplishment of LTE population) 
owing to the efficient charge-exchange reaction with H atoms, and (2) these lower and 
upper terms are strongly coupled with collisions which equalize the b values of these two 
terms. 

It was further found after several test calculations that this is a rather robust con- 
clusion, which may not be affected by any uncertainty in the adopted collisional cross 
section. That is, although little is known and considerable ambiguities are involved in 
the collisional rates due to H I atoms for this 1-2 transition (Cf_ 2 ; for which we used 
the scaled value estimated from Cf_ 2 ), we obtained essentially the same results even if 
Cj i _ 2 was completely neglected; i.e., only the Cf_ 2 itself is sufficiently large to bring the 
b vales of terms 1 and 2 into the same value of unity. 

We also carried out a simulation for the limiting test case of completely neglecting 
the collisional interaction between terms 1 and 2 (Cf_ 2 = Cj a _ 2 = 0). In this case, an 
appreciable overpopulation in the upper term 2 was observed (the extent of its NLTE 
departure growing toward the shallower layer) as a result of the downward cascade via 
the UV transitions connecting to term 2, which makes the NLTE line strength somewhat 
smaller compared to the LTE value (i.e., a positive correction in terms of A defined in 
Sect. 3.1) because of the resulting inequality of Sl/B > 1 in the upper layer. Even in 
such an extreme case, however, the extent of the computed A value is Js 0.05 dex in most 
range of the atmospheric parameters and definitely insignificant in the practical sense. 
Accordingly, we can confidently conclude that the non-LTE effect can be safely neglected 
for the [O i] lines, for which LTE must be a valid approximation. 
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4.4. [O/Fe] vs. [Fe/H] relation 

4.4.1. Procedure of reanalysis 

Now we are ready to study the [O/Fe] vs. [Fe/H] relation of metal-poor stars by reana- 
lyzing the published equivalent-width data of O I 7771/7774/7775 and [O i] 6300/6363 
lines while applying the non-LTE corrections based on our calculations. 

We adopted the same T e ff, logg, [Fe/H], and £ values as those used in the literature, 
from which the data of equivalent widths were taken. Regarding the model atmospheres, 
we used Kurucz's (1993) grid of ATLAS9 models (£ = 2 km s _1 ) which were interpolated 
with respect to T off , logg, and [Fc/H] of a star. As in Sect. 3.2, Kurucz's WIDTH9 
program was invoked for determining the LTE abundance (logeQ TE ) while using the line 
data given in Table 1. Then, the [O/Fe] ratio is obtained as 

[O/Fe] ee (log eg" 3 + A - 8.93) - [Fc/H], (4) 

where A is the non-LTE abundance correction, which was evaluated from W\, T e g, log g, 
and £ based on the procedure described in Sect. 3.3 [i.e., application of Eq. (1) with Table 
3] for the case of O I 7771/7774/7775 lines, and A = was adopted for the [O i] 6300/6363 
lines (cf. Sect. 4.3). Note that the solar oxygen abundance was assumed to be 8.93 (in 
the usual scale of log en = 12) according to Anders & Grevesse (1989). We treat [O/Fe] 
derived from O I 7771-5 lines and [O i] 6300/6363 lines separately, which we hereinafter 
referred to as [O/Fe] 77 and [O/Fe] 63, respectively. In case that equivalent-width data are 
available for two or three O 1 triplet lines, or for both [O 1] lines, we calculated [O/Fe] 
for each line and adopted their simple mean to obtain [O/Fe] 77 or [O/Fe] 63- 

4.4.2. Results of [O/Fe] 

We tried to use as many W\ data as possible taken from various sources, though our 
literature survey is not complete. Consequently, ~ 300 stars (~ 200 for [O/Fe] 77 and 
~ 100 for [0/Fe]63) covering a metallicity range of —3 & [Fe/H] <, could be analyzed, 
which may be sufficient for statistically meaningful discussion. The adopted data of the 
atmospheric parameters (T e s, log 9, £) the equivalent widths for each star are given in 
Table 5 (available only in electronic form), where the results of the LTE abundances 
and the non-LTE corrections/ abundances are also presented. The resulting [O/Fe] ratios 
(i.e., [0/Fe]77 or [0/Fe]63) plotted with respect to [Fc/H] for each star are shown in 
Fig. 8a, and the corresponding non-LTE corrections applied to the O 1 7771-5 lines are 
also plotted in Fig. 8b. The mean ([O/Fe] 77) and ([0/Fe] 63 ), averaged in each [Fe/H] 
bin of 0.5 dex, are given in Table 6. The spline-fit curves for the runs of ([0/Fe]77) and 
([O/Fe] 63) with [Fe/H] are also depicted in Fig. 8a. Our conclusions read from Fig. 8 and 
Table 6 are as follows. 
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(1) For mildly metal-deficient disk stars (-1 £ [Fc/H] & 0), our [0/Fe] 77 and [0/Fe] 6 3 
show a quite similar behavior; i.e., a gradual increase from [0/Fe]~ ([Fe/H] ~ 0) to 
[0/Fc]~ +0.5 ([Fe/H] ~ —1). This means that our non-LTE corrections are successful 
to achieve an agreement between the oxygen abundances derived from [O i] 6300/6363 
and O I 7771-5 for those population I stars, which is consistent with the conclusion of 
Takeda et al. (1998). 

(2) However, a discrepancy begins to appear when we go toward a lower metallicity 
further down from [Fc/H] ~ — 1. That is, [0/Fc]77 still increases with a slight change 

(i.e., gentler) of the slope from [0/Fe]~ +0.5 ([Fe/H] 1) to [0/Fc]~ +1 ([Fe/H] 

~ —3), while [0/Fe]@3 maintains nearly the same value at [0/Fc]~ +0.5 on the average 
over the metallicity range of —2.5 [Fe/H] ^ — 1 showing a quasi-flat plateau. This kind 
of [0/Fc]77 vs. [0/Fe]@3 discordance for those conspicuously metal-poor halo stars is just 
what has been occasionally reported so far (cf. Sect. 1). Our analysis thus suggests that 
this discrepancy in halo stars (amounting to ~ 0.2-0.3 dex at —2.5 [Fc/H] ^ —1.5) 
does exist, without being successfully removed by our non-LTE corrections, unlike the 
conclusion of Carretta et al. (2000). We will discuss this problem more in detail in the 
next section. 

4.5. O I vs. [O I] discrepancy 
4.5.1. Discordance in halo stars 

We have thus confirmed the existence of discrepancy between [0/Fe]77 and [0/Fc]63 in 
the quite metal-poor regime (—3 ^ [Fe/H] <, —1), while a good agreement was obtained 
in disk stars (-1 £ [Fe/H] £ 0). 10 

This can not be successfully removed by the non-LTE effect though they surely act 
in the direction of mitigating the discordance. According to our opinion (cf. Sect. 4.2), 
the O I vs. [O i] consistency obtained by Carretta et al. (2000) may be nothing but 
a fortuitous coincidence caused by their questionable non-LTE corrections for the O I 
7771-5 triplet (i.e., systematically overcorrected toward a lower metallicity). 

It is also unlikely that systematic errors in the T e ff scale (affecting mainly on O i) 
is responsible for the cause of this discrepancy as has been occasionally suspected (e.g., 
King 1993), since our analysis in Sect. 4.4 is based on the mixed data of equivalent widths 
and atmospheric parameters taken from a variety of data sources. 

Hence, we would consider that this discordance in halo stars is more or less inevitable 
as far as the conventional abundance determination method is concerned (i.e., "stan- 

Rather exceptionally, the recent data of Nissen et al. (2002) suggest more or less consistent 
abundances for both O I and [O i] (cf. double squares and circled pluses in Fig. 8) even for the 
very metal-poor stars as they confirmed in their ID analysis. We will return to this point in 
Sect. 4.5.2. 
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dard" non-LTE calculations as well as model atmospheres); that is, a challcngingly novel 
modelling would be required to bring them into agreement. 11 

4.5.2. Search for a trend in terms of T c g- and logg 

Then, the next question of interest is naturally "which solution is correct ?" , or in other 
words "which yields the wrong abundances, O I 7771-5 or [O i] 6300/6363 ?" . While this 
is a difficult question, it may be worthwhile as a first step to study the nature of the 
discordance; i.e., its dependence on T e g and log g. 

Let us invoke the following guiding principle: If we assume that the origin of the error 
is concerned with some stellar physical property which is not appropriately modelled, 
it is reasonable to expect that abundance discrepancies shown by the erroneous side 
(whichever O I or [O i]) may show some kind of systematic tendency in terms of these 
stellar parameters, since any stellar physical environment must undergo a distinct change 
if these parameters vary substantially. And, in an ideal case, the discordance might 
asymptotically disappear if we go to cither limit of the parameter variation. 

Motivated by this thought, we paid attention to [0/Fe]77 — ([0/Fe]63) (the discrepancy 
of each star's [0/Fe]77 relative to the mean ([0/Fc] 6 3}) and [0/Fc] 6 3 - ([0/Fe]77) (the 
inverse of the above) within two metallicity bins of —2.75 < [Fe/H] < —2.25 and —2.25 < 
[Fe/H] < —1.75, where a clear discordance of ([0/Fe]77) — ([0/Fc] 6 3) ~ 0.3dexis observed 
(cf. Table 6). The dependence of former quantity for each stars is plotted in Fig. 9 with 
respect to T c $ (a) and log g (b), while that of the latter is similarly shown in Fig. 10. 

In Fig. 9, what we observe is mostly a large scatter around the mean of <~ +0.3, which 
is presumably caused by the strong T e g-sensitivity of these high-excitation O I 7771-5 
lines. However, we can recognize in Fig 9a a marginal trend that [0/Fe]77 — ([0/Fe] 63 ) 
tends to decrease toward higher T e g (especially at T e g > 6000 K). Also, if we pay attention 
to the data for —2.75 < [Fe/H] < —2.25 in Fig. 9b (open circles), a slight log g-dependence 
appears to exist in the sense that higher-gravity stars tend to show less extent of the 
discordance. Similarly, Fig. 10 also suggests a weak trend for this group of very metal-poor 
stars (open circles); i.e., |[0/Fe] 6 3 — ([0/Fe]77)| appears to decrease with an increase in T e g 
as well as in log g toward the high-T ef j /high-log g limit. Consequently, a weak systematic 
trend may exist, especially for the very metal-deficient case of —2.75 < [Fe/H] < —2.25, 
in the sense that the discrepancy tends to be lessened for an increase in T c s as well as 
in log g. 

11 More realistic line- formation calculation based on 3D inhomogeneous model atmospheres 
might be a candidate. However, to our current knowledge, this effect does not improve the 
situation; rather, it even increases the discrepancy (i.e., the O I abundance is barely affected, 
while the [O i] abundance is somewhat lowered (cf. footnote 1). 
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However, these suspected T c ff- and log (/-dependences should not be regarded as being 
independent. We should recall that the very metal-poor stars ([Fe/H] —2) analyzed 
so far are mainly giants or subgiants. Namely, since T e g and logg show a close positive 
correlation in metal-poor giants reflecting the position on the HR diagram (see, e.g., Fig. 
3 of Pilachowski et al. 1996), the rather similar behavior of the O I — [O i] discrepancy 
in terms of the changes in both T c ff and logg may be naturally understood. 

Consequently, what we can learn from Figs. 9 and 10 is, modestly speaking, that 
analyzing higher T e s/\ogg subgiants (and maybe also dwarfs) has a less probability of 
suffering a large O I — [O i] discrepancy as compared to the study of lower T e ff / log g 
giants. The recent results of Nissen et al. (2001) (cf. footnote 10) may be regarded as 
supporting this speculation, since their data are based on stars with \ogg > 3. 

Unfortunately, however, Figs. 9 and 10 do not allow us any definitive arguments 
regarding the cause of the discrepancy or which of the O I or [O i] is problematic. It is 
apparent from these figures that the [O i] data of very metal-poor subgiants and dwarfs 
with higher T ff/ logg (and inversely the O I data for very metal-poor giants with lower 
T e ff/ log g) are quite insufficient, reflecting the preference of [O i] lines for giants so far 
(because their strengths increase with a lowering of the gravity). Hence, in order for a 
reliable statistical discussion, much more observations of [O i] lines for very metal-poor 
dwarfs or subgiants of higher gravity (along with O I lines for very metal-poor giants of 
lower gravity) are desirably awaited, which would surely improve this situation toward 
clarifying the origin of the discrepancy. 

5. Conclusion 

For the purpose of investigating the line formation of O I 7771-5 lines of IR triplet and 
[O i] 6300/6363 lines especially concerning the consistency of the oxygen abundances 
derived from these two in metal-poor disk/halo stars, we carried out non-LTE calculations 
with a new atomic model of oxygen on an extensive grid of model atmospheres (4500 K 
< T cff < 6500 K, 1 < log 5 < 5, and -3 < [Fc/H] < 0). 

In order to make the resulting non-LTE abundance corrections of O I 7771-5 lines 
more useful for practical applications, we derived an analytical formula with appropriate 
coefficients computed and tabulated in a grid of atmospheric parameters, by which the 
non-LTE correction is easily evaluated for any given set of W\, T s, logg and £. 

We discussed our non-LTE corrections for the O I triplet lines while comparing them 
with those derived from other studies published so far. We could not reproduce the 
trend of the non-LTE correction derived by Carretta et al. (2000), the systematic [Fe/H]- 
dependence of which is essentially the reason for their success of accomplishing the O I 
vs. [O i] consistency. 
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In contrast, we have shown that the non-LTE effect of O I 7771-5 lines is well described 
by the classical two-level-atom model and the metallicity can not be an essential factor. 
As a matter of fact, the formation of this O I triplet was found to be quite simple in the 
sense that only the transition-related quantities are important and details of the atomic 
model are not very significant, supporting the arguments of Kiselman (2001) and Nissen 
et al. (2002) who also arrived at the same conclusion. 

Regarding the forbidden [O i] lines, we arrived at a robust conclusion that the non- 
LTE effect for the [O i] 6300/6363 lines can be safely negligible and thus LTE is a valid 
approximation. 

We reanalyzed the various published equivalent-width data of O I 7771-5 and [O i] 
6300/6363 lines (while the non-LTE effect is taken into account for the former based 
on our calculations), in order to study the behavior of [O/Fe] vs. [Fe/H] relation over a 
wide range of metallicity (—3 & [Fe/H] ^ 0). While [0/Fe] 77 and [0/Fe] 6 3 were found 
to agree with each other for disk stars (—1 5s [Fe/H] 5s 0), the existence of a systematic 
discrepancy was confirmed for halo stars (—3 & [Fe/H] ^ —1) amounting to [0/Fe]77— 
[0/Fe] 63 ~ 0.3 dex at -3 <S [Fe/H] £ -2. 

An inspection of the dependence of discrepancy upon T e g and log g suggested that the 
extent of the discrepancy tends to be comparatively lessened for higher T c ff/log<7 stars, 
indicating the preference of dwarf (or subgiant) stars for studying the oxygen abundance 
of mctal-dcficicnt stars. 

Appendix A: On the effect and treatment of H i collision 

The inelastic collision rate due to neutral hydrogen atoms is regarded to be one of the 
important ambiguous factors in non-LTE calculations on late-type stellar atmospheres. 
A rough classical formula derived by Steenbock and Holweger (1984) is available for this 
(see Takeda 1991 for more detailed description on the standard classical recipe adopted). 
However, because of its nature of only an order-of-magnitude accuracy, a correction factor 
is often introduced, by which this classical rate is to be multiplied. We express this factor 
as 10 h by using the logarithmic correction (h). Negative h leads to the suppression of 
the H I collisional rates compared to the classical case, which increases the extent of the 
NLTE abundance correction (|A|); and vice versa. This situation is depicted in Fig. 11, 
where the /i-dependence of A for the O I 7774.18 line is shown, which was calculated for 
selected representative models. 

As mentioned in Sect. 2.1, we used h = (as was done by Takeda 1992, 1994) through- 
out this paper, which means the simple adoption of the classical treatment without any 
correction. This choice stems from the detailed investigation of Takeda (1995), where 
h = 0.0(±0.5) was concluded from the analysis of the solar flux spectrum based on (1) 
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the comparison of the solar oxygen abundance derived from O I 7771-5 lines with the 
average of other 15 weaker lines and (2) the profile-fitting of these near-IR triplet lines. 

It may be worthwhile, from another point of view, to check the consistency of this 
treatment on nearby solar- type stars, which have been recently analyzed by Takeda et 
al. (2001) to study the abundance characteristics of planet-harboring stars. 

Fig. 12 shows the comparison of the non-LTE abundance of O I 7771-5 triplet (with 
the NLTE correction obtained by using the formula described in this study) with the 
LTE abundance of O I 6158.2 (a) and [O i] 6300.3 (b), which were derived from either 
equivalent-widths or spectrum-synthesis technique (with the same line data as given in 
Table 1) as described in Takeda et al. (2001). Note that LTE can be safely applied to 
these latter two lines (cf. Table 7 and Sect. 4.3). 

We can see from Fig. 12a that loge^gg and log £7771—5 arc in satisfactory agreement 
with each other, which confirms that our choice of h = is surely reasonable at least 
within the framework of our non-LTE calculations. 

Meanwhile, for the case of 6300 vs. 7771-5 comparison shown in Fig. 12b, the scatter 
is comparatively large, which makes this figure less useful for the purpose of checking the 
adequacy of h. Presumably, this appreciable spread may be due to the large difference 
in the T e g -sensitivity of these two lines (in contrast to the case of 7771-5 and 6158 
which have almost the same T e fj-dependence). A close inspection of Fig. 12b suggests 
that the 6300 abundance tends to be somewhat larger especially for the case of higher 
oxygen abundance (log e ^ 9.0), for which two reasons may be suspected; i.e., the possible 
inaccuracy in the adopted gf value and/or the effect of blending. That is, according to 
the recent study of the solar [O 1] 6300 line profile by Allendc Prieto et al. (2001), the 
Ni 1 line at 6300.34 A makes an appreciable contribution to this feature. Our neglect of 
this effect may have cause an overestimation. In addition, their paper suggests that the 
most recently calculated gf value is logg/([0 1] 6300) = —9.72, which is by +0.1 dex 
larger than that adopted in this study (cf. Table 1). Hence, our solution of logeg 300 may 
have been somewhat overestimated also from this point of view. 

We finally remark that, even if our log eg 300 values have been overestimated either 
due to the use of underestimated log gf or due to the blending effect of the Ni line, our 
conclusion of 6300 vs. 7771-5 abundance discrepancy (cf. Sect. 4.4.2) can not be changed 
at all (i.e., because the discordance would become even larger). 

Appendix B: Interpretation of the triplet-line formation with the 
two-level-atom model 

It was shown in Sect. 4.2 that the formation of O 1 7771-5 lines can be well described 
by the "two-level" atomic model. Based on this simple line source function [cf. Eq. (2)] 
with the classical Milne-Eddington model (i.e., depth- independent line-parameters with 
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the continuum source function linearly increasing with r), we try to explain the results 
of our non-LTE calculations, that the NLTE correction for the O I 7773 triplet (for given 
T e g, log g, and £) is nearly a monotonic function of W\ (irrespective of any reasonable 
change in the metallicity) . Though this is a very simple schematic model, it gives us a 
deep insight on the formation of this O I line. 

B.l. What we can learn from the basic characteristics 

It was Hummer (1968) who first showed that the extent of the dilution in the two-level- 
atom line source function defined by Eq. (2) is determined by (a) the photon-destruction 
probability defined as the ratio of the radiative to collisional (downward) transition prob- 
ability (e' ~ C21/A21; cf. equation [11-7] in Mihalas 1978) and (b) the relative impor- 
tance of the thermalizing continuum determined by the line-to-continuum opacity ratio 
rj = kl/k c (cf. Hummer 1968). 12 Namely: 

(i) For the same e' (i.e., atmospheres with the same temperature and density), the extent 
of 77 essentially controls the dilution of Sl is (i.e. the larger 77, the more enhanced NLTE 
effect). 

(ii) On the other hand, the extent of the dilution in St naturally determines the impor- 
tance of the NLTE effect (i.e., the extent of the NLTE correction) as well as the strength 
of the observed line (i.e., equivalent width) in this case of scattering line formation. 
We now understand from "fact (ii)" that a close relationship exists between the strength 
of the O 1 7771-5 triplet lines and the corresponding non-LTE corrections. Meanwhile, 
"fact (i)" simply tells that the essential key parameter is r\ only (for the same e'), which 
further suggests that the difference in the metallicity does not play any essential role 
in the present problem. 13 Namely, as far as the line-to-continuum opacity ratio is the 
same (along with the same T e g and log g), the same extent of 5l, the same equivalent 
width, and the same non-LTE correction are expected irrespective of the metallicity, thus 
yielding the same A vs. W\ correlation. 

B.2. Numerical simulation toward a qualitative explanation of Fig. 1 

It may be instructive to compare the results from a simple simulation based on such a 
"two- level" model with the actually obtained A vs. W\ relation shown in Fig. 1. Though 
the value of e' is strongly depth-dependent in real atmospheres, we found e' ~ 1 (dwarfs) 
and e' ~ 0.1 (supergiants) in the line- forming region from actual model calculations. We 
adopted the Milne-Eddington model assuming the depth-independent line-to-continuum 

12 Hummer (1968) used a quantity called [3 (= kc/kl; the continuum-to-line opacity ratio) for 
the same purpose, which is just equivalent to rf 1 . 

13 While the possibilities of "indirect" metallicity effects via the temperature structure or the 
collisional rates surely exist, they are insignificant in the present case, as mentioned in Sect. 4.2. 
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opacity ratio r\ and the linear form of the Planck function, B(t) = Bo + Bit = Bq(1 + 
Bi/Bot) (t: continuum optical depth). We fixed Bq = 1 and changed the gradient 
B\{= Bi/ Bq). Again the determination of this gradient at 7773 A is rather difficult since 
it is not linear in terms of r in actual cases, but we tentatively chose by inspection of 
the model atmospheres B\ = 2 (low T e ff case) and B\ = 1 (high T e ff case). Hence, the 
radiative transfer problems (two-level-atom, complete frequency redistribution, and the 
pure Gaussian line profile) investigated by Hummer (1968) were numerically solved by 
using the Rybicki's (1971) scheme for various values of r\ ranging from 10 -2 to 10 3 for 
each of the three combinations of (B\, e'); (2, 1), (2, 0.1), and (1, 1). And then the 
emergent flux spectra were computed from the resulting solution of the source function 
(NLTE case), along with the LTE case of 5*l = B (i.e., the limit of e' — > oo), to obtain 
the equivalent width (W). 

In Fig. 13 are shown the LTE and NLTE curves of growth [panel (a)], and the corre- 
sponding NLTE correction A log 77 (= log ^nlte — logr/LTE; i-e., difference of the abscissa 
for a given W) as a function of Wnlte [panel (b)]. We can see from Fig. 13b that the 
computed A log 77 vs. Wnlte curves well resemble (at least qualitatively) the A vs. W\ re- 
lations observed in Fig. 1, in the sense that they can be approximated with the functional 
form of Eq. (1). 

It is worthwhile paying attention to the effects of changing e' or B\ observed in Fig. 
13b. As expected, the extent of the NLTE correction increases with a decrease in e'. 
Meanwhile, that | A. log 77 1 is larger for a smaller gradient (B\) of the Planck function may 
be attributed to the decreased sensitivity of Wlte (to a change in 77) owing to the lessened 
gradient (i.e., Wlte is relatively more affected by the Planck function B, in contrast to 
the NLTE case which is mainly controlled by the radiation field) . It should be noted that 
these two qualitative behaviors concluded from our simple model reasonably explain the 
tendencies observed in Fig. 1. That is, the larger |A| for a lowered gravity (at a given 
W\) is reasonably explained by the e'-dependence described above. Meanwhile, regarding 
the tendency of larger |A| for a higher T e g (at a given W\), the above-mentioned effect 
of decreased B\ may be responsible for it. 
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Table 1. Data of neutral oxygen lines adopted for abundance calculations. 
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Note — All data are were taken from Kurucz & Bell's (1995) compilation as far as available. 
Followed by first four self-explanatory columns, damping parameters are given in the last three 
columns: 

Gammar is the radiation damping constant, log7 rat j. 

Gammae is the Stark damping width per electron density at 10 4 K, log(7 e /AT e ). 

Gammaw is the van der Waals damping width per hydrogen density at 10 4 K, log(7 w /JVH). 

* Computed as default values in the Kurucz's WIDTH program (cf. Leusin, Topil'skaya 1987). 
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Table 3. Coefficients (a, b) for O I 7771.94, 7774.17, and 7775.39. 
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Note — The coefficients (a, b) used with the analytic formula [Eq.(l)] for evaluating the non-LTE correction, corre- 
sponding to each of the atmospheric models (T c ff and log g values shown in columns 2 and 3) for three values of the 
microturbulent velocity (§ — 1, 2, and 3 km s™ 1 ). The "W range" gives the minimum and the maximum value of the 
equivalent- width data (in mA) based on which the coefficients were established (i.e., indication of the range of the W\ 
values to which the formula with these coefficients can be applied). 
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Table 4. Non-LTE corrections of O I 7771.94 line for the same equivalent width but with 
different metallicities. 
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logr 
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Note — Shown are the NLTE corrections for the T c s = 6000 K and logg = 3.0 model at- 
mospheres of four different metallicities ([X] — [Fe/H]) of 0, —1, —2, and —3, but the oxygen 
abundances were carefully adjusted to yield the same O I 7771.94 equivalent widths W = 152 
mA (for [X] = and —1) and W = 20 mA (for [X] = —2 and —3) in order to enable the direct 
comparison with Gratton et al.'s (1999) results. r,? LTE is the residual flux at the line center and 
log r is the mean-depth of line formation [defined as J R\ log T5000 (ta = 2/3)dA/ / R x d\, where 
R\ (= l — r\) is the line-depth]; other symbols are self-explanatory. Note that our NLTE correc- 
tions (Aloge; column 11) are practically the same for the [X] = and —1 cases as well as for the 
[X] = —2 and —3 cases, in marked contrast with Gratton et al.'s (1999) results ( A log eGratton ; 
shown in column 12 for a comparison). 
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Table 6. [O/Fe] ratio averaged for each metallicity bin. 
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Note — 7177 and «63 are the numbers of stars used for calculating ([0/Fe]77) and {[0/Fe]63) at 
each metallicity bin, while 077 and 063 are the standard deviations. 
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Fig. 1. Non-LTE correction vs. equivalent-width relation calculated for O I 7774.17 (the 
middle line of the oxygen triplet) with the microturbulent velocity (£) of 2 kins" 1 , as 
constructed from (W / i NLTE , Aj) given in Table 2. Results corresponding to each of the 
models are discerned by differences in symbols, (a): high gravity cases (logg = 5.0, 4.0, 
and 3.0); (b): low gravity cases (log g = 3.0, 2.0, and 1.0). The curves drawn through the 
points show the analytical approximations [based on Eq. (1) with the coefficients given 
in Table 3], in which each line- type corresponds to the different surface gravity case: logg 
= 5 (dotted), 4 (dashed), 3 (solid), 2 (dashed), and 1 (dotted). 
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Fig. 2. Error of the analytical expression for the non-LTE correction (O I 7774.17, £ = 2 
km s _1 ), A = al0 bWx with the coefficients (a, b) given in Table 3, relative to the exactly 
calculated value (Table 2). Similarly to Fig. 1, the results are shown as functions of W\. 
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Fig. 3. Comparison of the non-LTE abundance corrections derived in previous two works 
of the author, with those calculated by our analytical formula [Eq. (1)] using appropriate 
(a, 6) coefficients evaluated by interpolating the values in Table 3. (a): O I 7771.94 in 
population I G-K giants (Takeda ct al. 1998); (b): O I 7771-5 in very metal-poor giants 
(Takeda et al. 2000). 
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Fig. 4. Comparison of our non-LTE abundance corrections for O I 7771-5 triplet result- 
ing from reanalysis of the published equivalent-width data by use of our approximate 
analytical formula with the original values determined by those authors, (a): Tomkin et 
al. (1992); (b): Mishenina et al. (2000); (c): Nissen et al. (2002). 



34 



Y. Takeda: Oxygen line formation in late-F through early-K disk/halo stars 




Fig. 5. (a): Comparison of our non-LTE abundance corrections for O I 7771-5 triplet 
resulting from reanalysis of the data of Tomkin et al. (1992; open circles) and Edvardsson 
et al. (1993; closed circles) (as was done by Carretta et al. 2000) with the corrections 
determined by Carretta et al. (2000). (b): Differences between Carretta et al.'s (2000) 
non-LTE corrections and ours, plotted as functions of metallicity. 
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Fig. 6. Non-LTE correction vs. equivalent-width relations for the O I 7771.94 line. Large 
symbols show the results calculated by Gratton et al. (1999), while small symbols indicate 
our results taken from Table 2. The model metallicity corresponding to each data point 
is denoted by the symbol type as explained in the figure. Our relevant (metallicity- 
independent) analytical relation [Eq. (1) coupled with Table 3] is depicted by the solid 
line, which was determined so as to fit the original data (small symbols) . (a) : T e f = 5000 
K, logg = 3.0; (b): T cf = 6000 K, logg = 3.0. 
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Fig. 7. Run of the physical quantities involved with the line formation of the O I 7771.94 
line, for the four demonstrative cases of the metallicity-dependence test given in Table 
4. The thin solid line and the thin dashed line indicate the depth-dependence of the line 
source function (in unit of the local Planck function) and the NLTE-to-LTE line-opacity, 
respectively (corresponding to the scales marked in the left-hand axis). The corresponding 
pure two- level-atom solution of S^/B [derived by solving the scattering problem with the 
line source function of Eq. (2) while using the actual depth-dependent e'(r), t](t), and 
B(t)] is also shown by the open circles, which is in good agreement with the full non-LTE 
solution presented by Eq. (3) (thin solid line). Meanwhile, the thick solid line and the 
thick dashed line show the run of the Planck function itself and that of the total (i.e., 
line plus continuum) source function at the line center, respectively (corresponding to the 
scales marked in the right-hand axis). The depth points corresponding to the continuum 
optical-depth of unity and the line-center optical depth of unity are indicated by two 
upward arrows. 
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Fig. 8. (a): [O/Fe] vs. [Fe/H] relation based on the reanalysis of published equivalent 
width data of O I 7771-5 lines (closed symbols, half-filled symbols, crosses) and [O i] 
6300/6363 lines (open symbols), where the non-LTE correction was taken into account 
only for the O I triplet lines. The solid and dashed lines indicate the mean [O/Fe] ten- 
dencies suggested from permitted O I lines and forbidden [O i] lines, respectively, (b): 
Applied non-LTE corrections for the O I 7771-5 lines plotted with respect to the metal- 
licity. Each of the symbols denote the sources of equivalent widths. O I 7771-5: filled 
squares — Carretta et al. (2000), St. Andrew's crosses — Edvardsson et al. (1993), filled 
inverse triangles — Abia & Rebolo (1989), Greek Crosses — Boesgaard & King (1993), 
half-filled triangles — Sncden et al. (1979), filled triangles — Mishenina et al. (2000), filled 
diamonds — Tomkin et al. (1992), filled circles — Cavallo et al. (1997), half-filled circles — 
Takeda et al. (2000), half-filled square — King (1994), double squares — Nisscn et al. 
(2002). [O i] 6300/6363: open squares — Carretta et al. (2000), open inverse triangles — 
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Fig. 9. Dependence of [0/Fe]777i_5 — ([0/Fc] 63 oo/6363) (the discrepancy of each star's 
[0/Fe]777i_ 5 relative to the mean ([0/Fe] 630 o/6363}) u P on ( a ) T cS and (b) logg. Open 
symbols are for stars of —2.75 < [Fe/H] < —2.25, while filled symbols correspond to those 
of -2.25 < [Fc/H] < -1.75. 
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Fig. 10. Dependence of [O/Fe] 6300/636 3 - ([0/Fe] 77 7i_ 5 ) upon (a) T off and (b) logg. 
Otherwise, the same as in Fig. 9. 
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Fig. 11. Dependence of the non-LTE correction for the O I 7774.18 line upon the loga- 
rithmic correction (h) applied to the classical treatment of H I collision rates. Left panel: 
([O/Fc] = 0.0, £ = 2 km s _1 ) calculations for representative [Fe/H] = models; right 
panel: ([O/Fc] = +0.5, £ = 2 km s _1 ) calculations for representative [Fe/H] = —2 models. 
Code "taagfo&mc" denotes the model with T e fj(K) = aa x 100, log g (cm s~ 2 ) = 66/10, 
and [Fe/H] (metallic! ty) = — c. 
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Fig. 12. Comparison of the oxygen abundances derived from the O I 7771-5 lines (with 
the NLTE corrections presented in this study) with those determined from other oxygen 
lines for nearby solar-type stars analyzed by Takeda et al. (2001): (a) O I 6158.2 line 
(comprising three components), (b) [O i] 6300 forbidden line. 
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Fig. 13. Results of the calculations on the simple "two-level atom" model for simulating 
the formation of the O I 7771-5 triplet. Shown are the results computed for three different 
combinations of B\ (gradient of the Planck function) and e' (photon destruction prob- 
ability) as indicated in the figure, (a) Theoretical curves of growth computed for LTE 
(without scattering) and NLTE (with scattering), (b) Run of the non-LTE correction 
(= log ?7nlte — logr/LTE) with the NLTE equivalent width. 



